Universal behaviour of entrainment due to coherent structures in turbulent shear flow 
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I suggest a solution to a persistent mystery in the physics of turbulent shear flows: cumulus clouds 
rise to towering heights, practically without entraining the ambient medium, while apparently similar 
turbulent jets in general lose their identity within a small distance through entrainment and mixing. 
From dynamical systems computations on a model chaotic vortical flow, I show that entrainment 
and mixing due to coherent structures depend sensitively on the relative speeds of different portions 
of the flow. A small change in these speeds, effected for example by heating, drastically alters the 
sizes of the KAM tori and the chaotic mixing region. The entrainment rate and, hence, the lifetime 
of a turbulent shear flow, shows a universal, non-monotone dependence on the heating. 



AZ 




FIG. 1. Schematic sectional view of a turbulent shear flow 
(a jet). Rotating the figure about the Z axis, the external 
coherent structures can be seen to be vortex rings. The small 
scale structures at the edge of the shear flow nibble at the 
ambient but do not entrain significantly. 

Turbulent shear flows [jjj in nature and in the labora- 
tory are anisotropic and inhomogeneous on large enough 
length scales |||). Among their most important prop- 
erties is entrainment Q, viz., the advection of the am- 
bient non-turbulent medium into the body of the turbu- 
lent flow, and the consequent expansion of the turbulent 
region. Coherent structures || |?j at the edges of the tur- 
bulent region are crucial to entrainment, but their role 
in turbulent transport in general is only now beginning 
to be understood |£j]. The entrainment process operates 
over three length scales ||. At the large scale, the co- 
herent structures 'engulf parcels of ambient fluid. The 
engulfed parcels then 'mingle' with original shear-flow 
fluid by breaking- up/stretching into smaller scales. Sub- 
sequently, the fluids mix at the molecular level when their 
characteristic sizes are comparable to molecular diffusion 
scales. 

The classic example of a turbulent shear flow is a jet 
(see Fig. |l|), which consists of a turbulent flow injected 



into a quiescent ambient fluid. Entrainment causes a jet 
to grow linearly in width, at a rate proportional to its 
axial velocity (Taylor's hypothesis pOfl ), as a function of 
position along its direction of propagation. The width of 
a plume - a hot jet - grows linearly as well, but faster 
than a jet Jllfl . It is therefore deeply puzzling that a cu- 
mulus cloud, formed by a rising, hot, moist plume, should 
entrain practically nothing, and rise to towering heights 
practically undiluted. Experiments and direct numerical 
simulations on cloud-like flows |l^,[f3| show, remarkably, 
that it is the local heating that is responsible for this 
shut-off of entrainment. This Letter shows that this in- 
triguing phenomenon can be understood by studying the 
Lagrangian trajectories of fluid particles in a dynamical- 
systems treatment of a model jet-like flow, and makes 
further important predictions. 

Our main results are as follows. We establish by quan- 
titative methods on a model that the entrainment and 
consequently the lifespan of a turbulent shear flow can 
be dramatically altered merely by tuning the relative ve- 
locities (buoyancy) between the coherent structures and 
the ambient. Surprisingly, a given buoyancy can either 
increase or decrease entrainment, depending on the spac- 
ing between the structures. Secondly, we show that the 
rate of engulfment is a universal function of the (appro- 
riately scaled) relative velocity. Although the inviscid 
model is incapable of estimating molecular mixing, we are 
able to make a good assessment of the mingling process 
using information about the chaos in the system. The re- 
sults lend themselves to verification by direct numerical 
simulations and by carefully designed experiments. 

I now motivate and construct the model that led to 
these results. The crucial step in entrainment is the 'en- 
gulfing' of ambient fluid by large-scale coherent struc- 
tures (Fig. |l|) which, as experiments as well as direct 
numerical simulations clearly show, are not destroyed by 
heating. It is thus vital to understand what heating does 
do. Such understanding comes not by solving for the 
flow behaviour in complete detail, but from an econom- 
ical model with only the essential features. To find such 
a model, note that the large-scale structures responsible 
for entrainment in a jet can be approximated (Fig. [l]) by 
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FIG. 2. Schematic of a Poincare section. The unstable 
manifold is shown by the solid line, the dotted line shows 
the stable manifold, and the crosses indicate a sample chaotic 
trajectory. The phase chosen is that at which the two rings, 
indicated here by diamonds, are concentric (at Z = 0). The 
frame of reference moves with the mean axial velocity of the 
rings. 

a stack of vortex rings. Vortical structures in the interior 
of the jet do not participate in the entrainment process 
p3| . The mixing features of such a flow are largely con- 
tained in an even simpler flow, that of a single pair of 
vortex rings. 

We begin with identical coaxial vortex rings of radius 
R = 1 and an axial separation AZ. Each ring moves 
forward with the sum of its self- induced velocity U s , the 
induced velocity Ui due to the other ring and the pre- 
scribed buoyancy velocity Ub- The rings are taken to 
be "heated/cooled" if they are buoyant relative to the 
surrounding fluid, instantaneous acceleration is neglected 
and the Boussinesq approximation used. The fluid mo- 
tion then proceeds as follows Jl4||l5| : for each vortex con- 
figuration, fluid lagrangian trajectories, obeying Hamil- 
tonian dynamics, may be obtained from the Biot-Savart 
law. The induced velocities from the vortices lead (a) to 
leap-frogging of the vortex rings, (b) to portions of the 
neighboring flow being trapped forever within KAM tori, 
and (c) to complex, often chaotic trajectories elsewhere, 
both within the jet and extending into the ambient. The 
dynamics in the last region is responsible for entrainment 
and subsequent mingling of ambient fluid with original jet 
fluid |]l6| , |l7|. These features may be seen on the Poincare 
map (Fig. |2|) vnr — Z space, where r and Z are the radial 
and axial coordinates respectively. The engulfmcnt may 
be obtained quantitatively |T(| from the tangling of man- 
ifolds in this map. The flow being axisymmetric, regions 
that appear as lobes in Fig. [| are funnel-shaped struc- 
tures of equal volume. During each half-cycle of leapfrog, 
the volume of fluid contained in one funnel (denoted by 
F in the figure) is entrained into the turbulent flow from 



1 







-3 




r 



FIG. 3. The unstable manifold for small ring spacing 
{AZ = 0.22). The stable manifold for U b = 4 as t -» oo 
is shown by the dotted line. (The self-induced ring velocity 
Us ~ 4.5, so U b /U s ~ 1). 
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FIG. 4. The unstable manifold for a large ring spacing 
(AZ = 1.00). 

the ambient, while an equivalent volume (D) is detrained 
from within the mixing region and dumped into the am- 
bient . The quantity engulfed in a fraction / of a cycle 
is equal to 2fF. The toroidal volume and the circulation 
of the rings are held respectively at 0.2ir 2 and 47r, fixing 
the basic velocity and time scales. 

In the absence of buoyancy, the present predictions 
agree with the results of p8| , corresponding to a small 
AZ, for a leap-frogging vortex pair. In addition, we find 
a maximum in the engulfment rate (proportional to the 
Melnikov function fl9|| ) as a function of the maximum 
separation AZ (inversely related to the frequency of leap- 
frog), which agrees with [^fj| . Our main results emerge 
when 'buoyancy' is switched on: for small ring spacing 
(AZ = 0.22), the entrainment is greatly enhanced (Fig. 
0) when the direction of buoyancy is such as to speed 
up the rings. In contrast, when the rings are slower, the 
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entrainment is killed off almost completely. Surprisingly, 
as AZ increases, the role of Ub is progressively reversed, 
until, at AZ = 1.0 (Fig. ||), the engulfment is inversely 
related to Ub- The extent of chaos within the region of 
influence of the vortical structures provides an estimate 
of the 'mingling' between original and ambient fluid. It 
is seen in Fig. || that for AZ = 0.22, when Ub — —2, a 
major portion of the region of interest consists of a se- 
ries of imbedded KAM tori, which make any degree of 
mingling unlikely. At Ub = — 1 (not shown) slow mix- 
ing takes place, while the chaotic nature of the flow at 
Ub = indicates that any entering fluid will follow a 
space-filling trajectory, undergo stretching into smaller 
scales and mingle with the original fluid, resulting in a 
dilution of the latter. Fig. |^ shows the dependence of F 
on the buoyancy for different values of ring spacing. For 
the highest spacing, due to the extreme tangling of the 
manifolds, the computations were reliable only for the 
range of Ub shown. When AZ is about half the radius, 
the engulfment decreases significantly whether the fluid 
is heated or cooled. For each ring spacing the engulfment 
displays a maximum at some value of the buoyancy ve- 
locity Ub- Table 1 shows characteristic axial velocities 
for different AZ. Here Ui is the maximum axial velocity, 
scaled by U s , induced by one of the rings on the other 
during a given cycle, (Ui would also be the typical ve- 
locity induced by one of the rings on a non- vortical fluid 
element at the same distance.) [C/f,] opt is the buoyancy 
velocity required to maximize engulfment at a given AZ. 
The last column in the table is independent of AZ, which 
suggests that we define a velocity ratio U r : 

Ur = (1) 

Ui + U b + cU s w 

whose numerator is a measure of the velocity of the in- 
duced motion in the surrounding fluid, and denominator 



represents the total ring speed. The constant c (1.55 
here) depends on the ratio of the ring thickness to the 
toroidal radius, as well as on the particular measures of 
the self and mutually induced velocities used. Since the 
period r of a leapfrog is vastly different for different ring 
spacings, the rate of engulfment E = F/t, would be a 
more appropriate measure for comparison. This quan- 
tity is plotted in Fig. as a function of the velocity ratio 
U r . 
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0.22 
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1.450 
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4.556 


1.56 
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1.236 


-0.327 


4.591 


1.56 


0.80 


1.109 


-0.433 


4.620 


1.54 



The rate of entrainment is seen to be a universal func- 
tion of the velocity ratio, i.e., the parameter which de- 
termines the entrainment rate is the ratio of the motion 
induced in the surroundings to the speed of the coherent 
structures. It is observed that the entrainment when the 
rings are much faster than the surroundings (extreme left 
of Fig. levels off at a small value. This feature will be 
investigated in detail elsewhere. 

In the simulations of Jl3| , the structure spacing is re- 
duced upon heating from about AZ — 1 to AZ = 0.3. 
The ratio U r in the present scales is about 0.45 in the 
pre-heating zone, where we expect a significant level of 
entrainment; and about 0.7 in the heating zone, where 
entrainment should drop almost to zero; which concur 
with what the simulations show. In a lab jet (or in a 
cloud), the entire turbulent flow would be heated with 



3 



6* 



0.5 



4 



V 



0.4 



\ 



AZ 
• 0.22 
<0.40 
□ 0.50 
a 0.65 
*0.80 



\ 

\ 



0.5 



0.6 



0.7 



U 



FIG. 7. Dependence of the entrainment rate on the velocity 
ratio. The fit -E/-B ma x = exp {— [(U r — £A. opt )/cr] 2 } is shown 
by the dashed line. The fitting parameters have the numerical 
values E max = 4.3 and a 2 — 0.09 here, which will change with 
the circulation prescribed. 



the system frequency. An increase in the engulfment of 
external fluid due to buoyancy is shown to be accompa- 
nied by an enlargement of the chaotic region, resulting 
in increased mingling. In closing, note that we have ne- 
glected viscous effects as well as the complicated manner 
in which heating/cooling is bound to affect a turbulent 
shear flow, by altering the geometry, density and orien- 
tation of the coherent structures, which would in turn 
affect entrainment. Including these would be natural di- 
rections for future research. In order to test the present 
predictions and to strengthen our understanding the ef- 
fect of heat, systematic experiments and direct numerical 
simulations are called for. The ideal experiment would be 
on a laboratory-generated isolated pair of vortex rings, 
whose circulations could be estimated and which could 
be heated or cooled over some axial distance. 

I thank A. Leonard, S. Wiggins, R. Narasimha and K. 
R. Sreenivas for discussions and useful suggestions. This 
work is funded by the AR&DB, India. 



respect to the ambient, not the structures alone. More- 
over, the rings would be of unequal sizes and strengths. 
We have examined these situations: when the entire in- 
ternal fluid, or a selected portion of it, is given the same 
buoyancy as the rings, the results are very close to those 
shown above, demonstrating that the important param- 
eter is the relative velocity between the structures and 
the external flow, and not the state of buoyancy of the 
interior fluid. Inequalities in the rings result in one ring 
becoming dominant over the other, but have very little 
qualitative effect on the results. 

In summary, it is known from observations that the 
entrainment in heated plumes is higher than in unheated 
jets, while that in (heated) clouds is much lower. This 
is the first computation to our knowledge that addresses 
this apparent contradiction, obtains entrainment quan- 
titatively albeit for a model heated flow, and shows 
that an imposed relative buoyancy can either increase 
or decrease the entrainment. For example in Fig. |t], 
when AZ = 0.5, the maximum entrainment occurs when 
Ub = 0, speeding-up or slowing-down the rings would 
both decrease entrainment. For small values of struc- 
ture spacing, the entrainment is directly related to the 
buoyancy (with the sign as defined here), while while at 
large spacing they are inversely related. The data col- 
lapse seen in Fig. [?] could not have been anticipated 
from intuition, and much less from an examination of 
the tangling of the manifolds. The funnel shapes and the 
convolutions they undergo are often qualitatively differ- 
ent in each case (they sometimes even have two or three 
sub- lobes per half-cycle), and so is the volume they en- 
train, but the rate of entrainment obeys a very simple 
relationship. Implicit in the data collapse is the remark- 
able observation that the maximum rate of engulfment 
(for a given circulation and geometry) is independent of 
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